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Abstract 

We study a static, spherically symmetric system of (2j + 1) mas- 
sive Dirac particles, each having angular momentum j, j — 1,2,..., 
in a classical gravitational and SU{2) Yang-Mills field. We show that 
for any black hole solution of the associated Einstein-Dirac- Yang/Mills 
equations, the spinors must vanish identically outside of the event hori- 
zon. 
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1 Introduction 



Recently the Einstein-Dirac- Yang/Mills (EDYM) equations were studied for 
a static, spherically symmetric system of a Dirac particle interacting with 
both a gravitational field and an SU{2) Yang- Mills field [1, 2]. In these 
papers, the Dirac particle had no angular momentum, and we could make 
a consistent ansatz for the Dirac wave function involving two real spinor 
functions. In the present paper, we allow the Dirac particles to have non- 
zero angular momentum j, j = 1,2,.... Similar to [3], we can build up a 
spherically symmetric system out of {2j + 1) such Dirac particles. In this 
case however, a reduction to real 2-spinors is no longer possible, but we can 
obtain a consistent ansatz involving four real spinor functions. 

We show that the only black hole solutions of our 4-spinor EDYM equa- 
tions are those for which the spinors vanish identically outside the black hole; 
thus these EDYM equations admit only the Bartnik-McKinnon (BM) black 
hole solutions of the SU{2) Einstein- Yang/Mills equations [4, 5]. This result 
extends our work in [2] to the case with angular momentum; it again means 
physically that the Dirac particles must either enter the black hole or escape 
to infinity. This generalization comes as a surprise because if one thinks of 
the classical limit, then classical point particles with angular momentum can 
"rotate around" the black hole on a stable orbit. Our result thus shows that 
the non-existence of black hole solutions is actually a quantum mechanical 
effect. A simple way of understanding the difference between classical and 
quantum mechanical particles is that for classical particles, the centrifugal 
barrier prevents the particles from falling into the black hole, whereas quan- 
tum mechanical particles can tunnel through this barrier. In our system, 
tunnelling alone does not explain the non-existence of black-hole solutions, 
because the Dirac particles are coupled to the classical fields; that is, they 
can influence the potential barrier. Our results are established by analyzing 
in detail the interaction between the matter fields and the gravitational field. 

In Section 2, we derive the static, spherically symmetric SU (2) EDYM 
equations with non-zero angular momentum of the Dirac particles. By as- 
suming the BM ansatz for the YM potential (the vanishing of the electric 
component), the resulting system consists of 4 first-order equations for the 
spinors, two first-order Einstein equations, and a second-order equation for 
the YM potential. This EDYM system is much more complicated than the 
system considered in [2], and in order to make possible a rigorous mathe- 
matical analysis of the equations, we often assume (as in [3]) a power ansatz 
for the metric functions and the YM potential. Our analysis combines both 
geometrical and analytic techniques. 
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2 Derivation of the EDYM Equations 

We begin with the separation of variables for the Dirac equation in a static, 
sphericahy symmetric EYM background. As in [1], we choose the hne ele- 
ment and the YM potential A in the form 

ds'^ = -i— dt'^ - - d-d^ - sin^ ^ d(p'^ (1) 

r(r)2 A{r) ^ ^ ' 

A = w{r)T^d-& + {cos-&T^ + w{r) sini?r^)d(^ (2) 

with two metric functions A, T, and the YM potential w. The Dirac operator 
was computed in [1, Section 2] to be 

G = iT^'dt + f [i^dr + J (vC4- 1) - ^ VI |J) 

+ i-i'^d^ + i-i^d^ +'^{w-l) (7f - fr') . (3) 

This Dirac operator acts on 8-component wave functions, which as in [1] 

we denote by ^' = (^""")a,«,a=i,2) where a are the two spin orientations, u 
corresponds to the upper and lower components of the Dirac spinor (usually 
called the "large" and "small" components, respectively), and a is the YM 
index. The Dirac equation is 

(G-m)* = 0, (4) 

where m is the rest mass of the Dirac particle, which we assume to be positive 
(m > 0). 

As explained in [1], the Dirac operator (3) commutes with the "total 
angular momentum operators" 

J = L + 5 + f , (5) 

where L is angular momentum, 5* the spin operator, and r the standard 
basis of su(2)ym- Thus the Dirac operator is invariant on the eigenspaces of 
total angular momentum, and we can separate out the angular dependence 
by restricting the Dirac operator to suitable eigenspaces of the operators 
J. Since (5) can be regarded as the addition of angular momentum and 
two spins \, the eigenvalues of J are integers. In [1], the Dirac equation 
was considered on the kernel of the operator J^; this leads to the two- 
component Dirac equation [1, (2. 23), (2. 24)]. Here we want to study the 
effect of angular momentum and shall thus concentrate on the eigenspaces of 
with eigenvalues j(j + 1), j = 1, 2, . . . . Since the eigenvalues of merely 
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describe the orientation of the wave function in space, it is furthermore 
sufficient to restrict attention to the eigenspace of Jz corresponding to the 
highest possible eigenvalue. Thus we shall consider the Dirac equation on 
the wave functions with 

J2*=j(j + 1)* and J,*=j^' (j = l,2,...). (6) 



Since (6) involves only angular operators, it is convenient to analyze these 
equations on spinors ip) on S"^ (where a and a are again the spin and 

YM indices, respectively). Let us first determine the dimension of the space 
spanned by the vectors satisfying (6). Using the well-known decomposition 
of two spins ^ into a singlet and a triplet, we choose a spinor basis ^st with 
s = 0, 1 and —s<t<s satisfying 

{S + ff ^st = sis + 1) ^st , (Sz + Tz) ^st = t^st- 

The spherical harmonics (^ifc)/>o, -i<fe</) on the other hand, are a basis of 
L^{S'^). Using the rules for the addition of angular momentum [6], the wave 
functions satisfying (6) must be linear combinations of the following vectors, 

Yjj^oo (7) 

^.-i.-i^ii (8) 

Yjj-i^ii , Yjj^.o (9) 

Yj+ij^io, yj+ij+i$i_i. (10) 

These vectors all satisfy the second equation in (6), but they are not nec- 
essarily cigcnfunctions of J^. We now use the fact that a vector ^' 7^ 
satisfying the equation J^^* = j^' is an cigcnstatc of with eigenvalue 
j[j + 1) if and only if it is in the kernel of the operator J"*" = Jx + iJy 
Thus the dimension of the eigenspace (6) coincides with the dimension of 
the kernel of J+, restricted to the space spanned by the vectors (7)-(10). A 
simple calculation shows that this dimension is four (for example, wc have 
J+ {Y^j j-i ^11) = Yjj = J+ {Yjj $10), and thus J+ applied to the 
vectors (9) has a one-dimensional kernel). 

We next construct a convenient basis for the angular functions satisfying 
(6). We denote the vector (7) by $o- It is uniquely characterized by the 
conditions 

L^^o = j(i + l)^o, Lz^o = j<^o 

(5 + f)$o = , ||$o||s2 = 1. 



We form the remaining three basis vectors by multiplying $0 with spheri- 
cally symmetric combinations of the spin and angular momentum operators. 
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namely 



$1 = 25'' $0 = -2r''$o 

$2 = - (SL) $0 = -- (tL) $0 
c c 

$3 = -5V5L)$o = --T''(fL)$o 
c c 



where 



Since the operators S'^' , t^', and (SL) commute with J, it is clear that the vec- 
tors $1,... ,$3 satisfy (6). Furthermore, using the standard commutation 
relations between the operators L, x, and S [6], we obtain the relations 

4 2 

(5t)2 = i - i S^f = - - - 5f 
^ ^ 4 2 16 2 

25'' $0 = -2r''$o = ^1 

25'' $1 = -2t''^i = $0 

25^ $2 = 2r^$2 = $3 
25'' $3 = 2r''$3 = ^2 

(5f)cI.o = $0 = -^$0 

(5t) $1 = 2Sfcrik S'' $0 = -2Sk (xS) Sk $0 

(5f)$2 = -{Sf)iSL)<^>o = -{LT)'^o--{SL){Sf)'^>o 
c c c 

= -- {SL) $0 + 1- (SL) $0 = 7 $2 
c zc 4 

(St) $3 = - (St) S''(Sr) $0 = - {SL) - - S'^ {St)^ $0 
c c c 

(SL)$o = 1*2 

(5L)$i = 2(SL)S'"$o 

= 2 Sj [Lj, S^] $0 + 2 {Sj, S^} Lj ^0 - 2S^ {SL) $0 

= -2i Sj ejki Xk Si^o - ^ $3 
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{SL) $2 
(SL) $3 



c c 
-ejki Xj ekim 'S'm $0 - 2 ^3 = -^1 - 2 



- {SL f $0 
c 



c V4 



\SL]^0 



$3 

2 ^ 



2 



2 " 



- $2 

2 



^ (SL) {SL) $0 = -5'' (SL) ^0 - cS'' $o 



and thus 



2(Sr-5V)r'' $o 

2(5f-5V")r''$i 

2(5f-5V")r'' $2 
2{St-SW')t'' $3 



2 

2 " 


. 



Using these relations, it is easy to verify that the vectors "i&O) • • • 1^3 are 
orthonormal on L^(S'^). We take for the wave function ^ the ansatz 



+ iP{r) ^ri^, ^) <2 + iSir) ^) K2) (11) 

with real functions a, /?, 7, and (5, where > is the energy of the Dirac 
particle. This is the simplest ansatz for which the Dirac equation (4) reduces 
to a consistent set of ODEs. Namely, we obtain the following system of ODEs 
for the four-component wave function $ := (a, /3, 7, 5), 



w 
r 

LoT-m 
c 

r 




-uT-m 
w 

r 








c 
r 



r 




ooT — m 



r 

-uT-m 



(12) 



Similar to [3], we consider the system of (2j + 1) Dirac wave functions 
obtained from (11) by applying the ladder operators J±. Substituting the 
Einstein and YM equations [1] and using the ansatz (1) and (2), we get the 
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following system of ODEs, 



2rA 



T' 



2A w' 



,/2 



(13) 



2 A w' 



/2 



4c 



4w; 



2m (a^ -P^ + Y-S^) + — (ad + ^7) + — a/3 

r r 



(14) 



-w(l-u;^) + e^rTap - ^r'^ A' w' + . (15) 



Here (13) and (14) are the Einstein equations, and (15) is the YM equation. 
Notice that the YM equation does not depend on 7 and S; moreover the lower 
two rows in the Dirac equation (12) are independent of w. This means that 
the Dirac particles couple to the YM field only via the spinor functions a and 
p. Indeed, a main difficulty here as compared to the two-spinor problem [2] 
will be to control the behavior of 7 and S. 

For later use, we also give the equations for the following composite 
functions, 



{Aw'y 
r {AT^y 



-w{l-w'') + e^rTaP + - r 



1 „2 {AT^y , 



2^2 



w 



2m (a^-/3^+7^-(5^) + — {a5 + (3-^) + — a(5 

r r 



(16) 



.(17) 



Also, it is quite remarkable and will be useful later that for cj = 0, the 
squared Dirac equation splits into separate equations for (0,7) and {P,6); 
namely from (12), 



VJdriVAdr^) 



m + 







( if)' 





c 





^ 




( 


nP 





wc 





\ 


+ 







-(7)' 





c 




+ 







7T 





wc 




c 















wc 

















V 


c 

f7 












[ 





wc 
7T 








/ 



$ . (18) 



1238 



ABSENCE OF BLACK HOLE SOLUTIONS 



3 Non-Existence Results 

As in [2], we consider the situation where r = p > is the event horizon of 
a black hole, i.e. A{p) = 0, and A{p) > ii r > p. We again make (cf. [1]) 
suitable assumptions on the regularity of the event horizon: 

(Ai) The volume clement -^1 det gij\ = \sm-d\ r'^A"^ is smooth and 
non-zero on the horizon; i.e. 

T-'^A-\ T"^ A eC\[p, oo)) . 
(An) The strength of the Yang- Mills field Fij is given by 

TV(F,,F^^) = — ^ + . 

We assume that it is bounded near the horizon; i.e. 

w and Aw''^ are bounded ioi p < r < p + s. (19) 

Furthermore, the spinors should be normalizable outside and away from the 
event horizon, i.e. 

OO rp 

|$P —= < OO for every £ > 0. (20) 

p+s V A 

Finally, we assume that the metric functions and the YM potential satisfy a 
power ansatz near the event horizon. More precisely, setting 

u = r — p , 

we assume the ansatz 

A{r) = Aqu' + o{u') (21) 
w-wq = wiu'^ + o{u'^) (22) 

with real coefficients Aq and wi, powers s,k > and wq = linir\p'w{r). 
Here and in what follows, 

f{u) = o{u'^) means that 36 > with limsup \u~^~^ /(■'^)| < oo . 

r\p 

Also, we shall always assume that the derivatives of a function in o{u'') have 
the natural decay properties; more precisely, 

fiu) = oiu") implies that /("^(^x) = o{u''-'') . 
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According to Ai, (21) yields that T also satisfies a power law, more precisely 

T(r) ~ u'i + o(n^5) . (23) 
Our main result is the following. 

Theorem 3.1. Under the above assumptions, the only black hole solutions 
of the EDYM equations (12)-(15) are either the Bartnik-McKinnon black 
hole solutions of the EYM equations, or 

4 1 
s = - and K = - . (24) 

In (24), the so-called exceptional case, the spinors behave near the horizon 
like 

{a/3){r) ~ u^+o{u^), < {jS){p) < 00. (25) 



Our method for the proof of this theorem is to assume a black hole 
solution with $ ^ 0, and to show that this implies (24) and (25). The proof, 
which is split up into several parts, is given in Sections 4-7. 

In Section 8, we will analyze the exceptional case. It is shown numerically 
that the ansatz (24), (25) does not yield global solutions of the EDYM equa- 
tions. Prom this we conclude that for all black hole solutions of our EDYM 
system, the Dirac spinors must vanish identically outside of the event hori- 
zon. 



4 Proof that u = 



Let us assume that there is a solution of the EDYM equations where the 
spinors are not identically zero, $ ^ 0. In this section we will show that 
then LJ must be zero. First we shall prove that the norm of the spinors |^>| 
is bounded from above and below near the event horizon. We distinguish 
between the two cases where A~2 is or is not integrable near the event 
horizon. 

Lemma 4.1. If A~2 is integrable near the event horizon r = p, then there 
are positive constant k and e such that 

7 < Mr)f < k, if p<r < p + e. (26) 
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Proof. Writing (12) as -s/A = M$, we have 

2 dr' ' 2 ^ ^ 

= - (a^ - /32) _ 2m (aP + -fS) + — (aj - m 
r r 

r r 
< ci |$p . (27) 

Here the constant ci is independent of r G (p, p+l], since w is bounded near 
the horizon according to assumption An. Since we are assuming that ^> ^ 
in r > p, the uniqueness theorem for solutions of ODEs yields that |$p > 
on {p,p + 1]. Then dividing (27) by ^\/]4|$p and integrating from ri to r2, 
/9 < ri < r2, we get 

|log|$(r2)P - log|$(ri)p| < 2ci / A~^{r)dr. 

J n 

Taking the limit ri \ p in this last inequality gives the desired result. ■ 



Lemma 4.2. If A 2 is not integrable near the event horizon r 
CO ^ 0, then there are positive constants k and e such that 



- < |$(r)|2 < k 



ifp<r<p + e. 



p and 



(28) 



Proof. Define the matrix J by 



w 



J 



ruT 


c 



w 

~ruT 
m 



\ ruT 



rcoT 




ruT 
m 







ruT 




1 



m 



and notice that, since r(r) — ^ 00 as r \ p, J is close to the identity matrix 
for r near p. If we let 

F{r) = <#(r), J(r) $(r)> , 
then a straightforward calculation yields that 

F' = <$(r), J'(r) $(r)> . 
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In a manner similar to that in [2], we ean prove that \J'\ is integrable near 
r = p, and as in [2], it follows that (28) holds. ■ 



Lemma 4.3. // $ ^ for r > p, then a; = 0. 



Proof. Assume that a; 7^ 0. We write the (AT^)' equation (17) as 



r {AT^ 
+ 



-4a; |$| 



4Aw'^ 



Aw 



4c 



2m (a^ -P^ + Y-S^) + — a(3+ — (aS + /?7) 



(29) 



According to hypothesis An , the left side of this equation is bounded near 
the event horizon. The Lemmas 4.1 and 4.2 together with An imply that 
the coefficients of T^, T^, and in this equation are all all bounded, and 
that the coefficient of is bounded away from zero near r = p. Assumption 
An implies that T{r) — > 00 as r \ p. Hence the right side of (29) diverges 
as r \ p. This is a contradiction. ■ 



5 Reduction to the Case a{p) = 0, (3{p) ^ 

Since a; = 0, the Dirac equation (12) reduces to 



/ 


w/r 


—m 


c/r 





\ 




—m 


—w/r 





—c/r 






c/r 








—m 




V 





—c/r 


— m 





/ 



(30) 



The following Lemma gives some global information on the behavior of the 
solutions to (30). 

Lemma 5.1. The function {oi(5 + 7(5) is strictly positive, decreasing, and 
tends to zero as r ^ oo. 

Proof. A straightforward calculation gives 

\/~A{af3 + -^5)' = -m|$p, 

so that (a/3 + ^5){r) is a strictly decreasing function, and thus has a (pos- 
sibly infinite) limit as r — oo. Since > 2 \af5 + ^5\, we see that the 
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normalization condition (20) holds only if this limit is zero. It follows that 
{aP + jd) is strictly positive. ■ 



Next we want to show that the spinors have a (possibly infinite) limit as 
r \ p. When A~2 is integrable near the event horizon, it is an immediate 
consequence of Lemma 4.1 that this limit exists and is even finite. 

Corollary 5.2. If A~2 is integrable near the horizon, then $ has a finite 
limit for r \ p. 

Proof. We can integrate (30) from ri to r2, p < ri < r2, 

n 1 

$(r2)-$(ri) = A 2{r) M{r) ^{r) dr . 

Lemma 4.1 yields that the right side converges as r\ \ p, and hence $ has 
a finite limit. ■ 

In the case when A~^ is not integrable near the horizon, we argue as follows. 
According to the power ansatz (22), the matrix in (30) has a finite limit on 
the horizon. Exactly as shown in [7, Section 5] using the stable manifold the- 
orem, there are fundamental solutions of the Dirac equation (that is, a basis 
of solutions of the ODE (30)) which behave near the event horizon exponen- 
tially like exp(Aj- / A~2)^ where \j G IR are the eigenvalues for r \ p of the 
matrix in (30) (notice that the \j are real since they are the eigenvalues of 
a symmetric matrix). Thus for any linear combination of these fundamental 
solutions, the spinor functions are monotone in a neighborhood of the event 
horizon, and hence as r \ /O, $ has a limit in IR U {±00}. We set 

$(p) = lim$(r), (q;/3)(p) = lim(a/3)(r) . 

r\p r\p 

Proposition 5.3. {aP){p) = 0. 

Proof. We consider the {Aw')' equation (16) with a; = 0, 

r' {Aw'y = -wil - w') + r (VAT) + '-^^^ Aw' . (31) 
Suppose that 

{a(3){p) > . (32) 
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From hypotheses Ai and and An, we se that the coefficient of a(5A~2 is 
positive near r = p, and the other terms on the right side of (31) are bounded. 
Thus we may write (31) in the form 

{AwJ = $(r) + 4== , (33) 

where # is bounded and ^' > near p. Thus we can find constants ^q, ^'q 
satisfying 

{Aw'Y > $0 + ^== , *o > 0, (34) 

for r near p. Then exactly as in [2, Section 3], it follows that the spinors 
must vanish in r > p. 

If on the other hand 

iaP)ip) < . (35) 
then (33) holds with ^'(r) < near p. Thus 

-{AwJ = -$(r) - 4= • (36) 

Setting w = —w, (36) becomes 

{Aw'y = -$(r) - ^^^'■^ 



where — ^'(r) > for r near p. Thus we see that (33) holds for w replaced 
by w. This again leads to a contradiction. ■ 



The next proposition rules out the case that both a and /3 vanish on the 
event horizon. 

Proposition 5.4. Either a{p) = 0, /3(p) ^ or a{p) ^ 0, /3(p) = 0. 
Proof. Suppose that 

aip) = = Pip) . (37) 

According to Lemma 5.1, 7 and 6 cannot both vanish on the event horizon. 
Using (30), we have for r near p, 

^ol = -7 + o(l) (38) 
r 

^/A(3' = --S + o(l) . (39) 
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HA 2 is not integrable near the event horizon, these equations show that 
7(p) and 5{p) are finite (otherwise multiplying (38) and (39) by A~2 and 
integrating would contradict (37)); if A~2 is integrable near p, Corollary 5.2 
shows that ^{p) and S{p) are again finite. 



From (17) with u; = we have 
r{AT'^y = 



A on An 

r r 



- ^ {Aw'"") . 



(40) 



Since the coefiicients of and are bounded, as is the left-hand side, we 
conclude that, since T{r) oo as r \ p, the coefficient of must vanish 
on the horizon. 



4w 



4c 



2m (a^ -f3^ + r-S ) + — a(3 + — (ad + 







r=p 



As a consequence, 7(/9)^ = S{p)^, and Lemma 5.1 yields that 

lip) = S{p) ^ . 
Furthermore from (38) and (39), for r near p, 

sgn Q;(r) = sgn 7(r) and sgn /3(r) = — sgn b{r) 



(41) 

(42) 
(43) 



From (42) and (43), we sec that for r near /), the spinors must lie in 
the shaded areas in one of the two configurations (I) or (II) in Figure 1. 
Now we claim that in either configuration (I) or (II), the shaded regions are 
invariant. For the proof, we consider the Dirac equation (30). One easily 
checks that the shaded regions in the a//?-plots arc invariant, provided that 
7 and 5 are as depicted in their shaded regions. Similarly, one verifies that 
the shaded regions in the 7/(5-plots are invariant, provided that a and /3 lie 
in the shaded regions. Moreover, Lemma 5.1 shows that the spinors cannot 
leave their regions simultaneously (i.e. for the same r). This proves the 
claim. 

Next we consider the situation for large v. In the limit r — oo, the 
matrix M in (30) goes over to the matrix S given by 

/ -m \ 
-m 

-m ■ 

V -m / 
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Figure 1: Invariant Regions for the Spinors 



In S, the non-zero 2x2 upper and lower triangular blocks, 

/ -m\ 

\ -m y ' 

have eigenvectors (1,1)* and (1,-1)* with corresponding eigenvalues —m 
and m, respectively. Since the system of ODEs 

splits into separate equations for (a, (3) and (7,5), we see that {a{r), f3{r)) 
must be a linear combination of e^'^^^^ ^ (1, 1)* and e'^*-''^ ^ (1, —1)*, where the 
functions c and d are close to m. Since the spinors are assumed to be nor- 
malizable (i.e. (20) holds), and are non-zero for r > p, it follows that for 
large r, the spinors are close to a constant multiple of e""^^'')'' (1,1)*, and 
thus for large r, sgn a{r) = sgn (3{r). Similarly, for large r, sgn j{r) = 
sgn 6{r). This is a contradiction to the shaded invariant regions of Figure 1. 



The two cases in Proposition 5.4 can be treated very similarly. Therefore 
we shall in what follows restrict attention to the first case. Furthermore, we 
know from Lemma 5.1 and Proposition 5.3 that ('yS){p) > 0. Using linearity 
of the Dirac equation, we can assume that both j{p) and 6{p) are positive. 
Hence the remaining problem is to consider the case where 



a{p) = , Pip) , 7(P), S{p) > . 



(44) 
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6 Proof that A 2 is Integrable Near the Event Hori- 
zon 



In this section we shall assume that A~ 2 is not integrable near the event hori- 
zon and deduce a contradiction. We work with the power ansatz (21), (22) 
and thus assume that s > 2. 

We first consider the case wq ^ 0. The first component of the squared 
Dirac equation (18) is 



VJdri^/Adra) 



rr? + 5 — + - 



« + (45) 



The square bracket is bounded according to An. Since a{p) = and 7(p) > 
0, our assumption wq ^ implies that the right side of (45) is bounded away 
from zero near the event horizon, i.e. there are constants 6, e with 

±VAdr{VAdra) > 6 ioi p < r < p + e, 

where "it" corresponds to the two cases wo > 1 and u;o < 1, respectively. 
We multiply this inequality by A~2 and integrate from ri to r2, p < ri < r2, 

r2 



r- '•2 1 

±\fAdra > 5 / A- 2 

^1 Jri 



The right side diverges as ri \ p, and thus limr\p ^/Adj-a = =Foo. Hence 
near the event horizon, ^dra > A" 2^ and integrating once again yields that 
lim^x^pa = ±00, in contradiction to a{p) = 0. 

Suppose now that u^o = 0. We first consider the A-equation (13), which 
since a; = becomes 



rA' = l-A-LJlz^.l^^„^^. ,46) 



Employing the power ansatz (21), (22) gives 

0(u'-^) = 1 + 0(u') - ^ + Oiu'^'') + 0(u^+2'^-2) . (47) 

e r 



Here and in what follows, 

f{u) = 0{u'^) means that lim u~'^ f{u) is finite and non-zero, 
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also we omit the expressions "o(u^)." The constant term in (47) must vanish, 
and thus e^p^ = 1. Using also that 0{u'^) is of higher order, (47) reduces to 

0{u'-^) = 0{u) + Oiu^") + 0{u'+^''-^) . (48) 

Suppose first that s > 2. Then (48) yields that k = ^. Substituting our 
power ansatz into the 74'u;'-equation (16) gives 

and thus a/? = 0{u~). Since P{p) ^ 0, we conclude that there are 
constants ci , 5 > with 

l+s 

\a\ < ci u~ ioT p < r < p + S. (49) 

From this one sees that the first summand on the right side of (45) is of 
higher order; more precisely, 

VAdriVAdra) = 0{u^) . 

Multiplying by A~2 and integrating twice, we conclude that 

a = 0{ui-') , 

and this contradicts (49). 

The final case to consider is u;o = ^'^d s = 2. Now the ^lu'-equation 
(16) gives 

0{u'^) = 0{u'') + e'^TaP 

and thus a = o{u). This gives a contradiction in (45) unless w — \fA = 
o(u), and we conclude that k = 1. Now consider the Dirac equation (30). 
Since w{p) = 0, the eigenvalues of the matrix in (30) on the horizon are 
A = ziz^/m^ + (? I p^. As a consequence, the fundamental solutions behave 
near the horizon ~ y±y/rn?+c?T^ _ boundary conditions (44) imply that 
a ~ y,fr\/^^^^r^l^ whereas /3, 7, 5 ~ ^^-aA^^+cV?^ and we conclude that 

(a(5 + /?7)(p) > 0. (50) 

Next we consider the ^T^-equation (17), which for a; = takes the form 
(40). It is convenient to introduce for the square bracket the short notation 

4.7/1 Af> 

[] = 2m(p? - P^ ^-i^ - S'^)^ — oiP + — {ad + p-f) . (51) 
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We define the matrix B by 



B 



/ m 


w/r 





c/r 


\ 


w/r 


—m 


c/r 










c/r 


m 







V c/r 








— m 


/ 



A short calculation shows that 

[] = 2<*, , 
and furthermore, using the Dirac equation (30), 

[]' = 2<^, B'^> = 4ap(^y - ^{aS + Pj) 



4w' ^ 4.W ^ 4.C , . ^ . 
ap ap K [ad + p^) . 



(52) 



Since {al3){p) = and {a5 + P^){p) > according to (50), 

— []'> C2 ioip<r<p + 6 
and a constant C2 > . Integrating on both sides shows that 

I [ ] I > C3 « iov p <r < p + 6 

with C3 > 0. As a consequence, the first summand in (40) diverges for r \ p, 
whereas the left side and the second summand on the right are bounded in 
this hmit. This is a contradiction. 

We conclude that A~2 must be integrable near the event horizon, and 
so s < 2. 



7 Proof of the Main Theorem 



In this section we shall analyze the EDYM equations with the power ansatz 
(21), (22) near the event horizon. We will derive restrictions for the powers 
s and k until only the exceptional case (24) of Theorem 3.1 remains. So 
far, wc know from Section 6 that s < 2. A simple lower bound follows from 
the A-equation (13) which for to = simplifies to (46). Namely in view 
of hypothesis An, the right-hand side of (46) is bounded, and thus s > 1. 
The case s = 1 is excluded just as in [2] by matching the spinors across 
the horizon and applying a radial flux argument. Thus it remains only to 
consider s in the range 

1 < s < 2 . (53) 
We begin by deriving a power expansion for a near the event horizon. 
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Lemma 7.1. Suppose that wq ^ or n ^ s/2. Then the function a behaves 
near the horizon as 

a = aou'^ + o{vF) , ao 7^ 0, (54) 

where the power a is either 

a = l-| (55) 

or 



2 — s if vjQ 

2 — s + min(/c, s/2) if wq = 0. 



(56) 



Proof. We set 



a = sup < p : limsup|u ^ a{r)\ < oo > < oo . (57) 
[ r\p J 

Suppose first that a < oo. Then for every ly < a there are constants c > 
and £ > with 

|a(r)| < cu^ for p < r < p + e. (58) 

We consider the first component of the squared Dirac equation (45) and 
write it in the form 

VAdr{\/^dra) = fa + g, (59) 

where / stands for the square bracket and g for the last summand in (45), 
respectively. Multiplying by A~2 and integrating gives 

^/A^ra{r) = [ A'^ (fa + g) + C 
J p 

with an integration constant C. We again multiply by A~^ and integrate. 
Since a(p) = 0, we obtain 

a{r) = j A-^{s)ds j A"^ {fa + g) + cf A'^ . (60) 

J p J p J p 

Note that the function /, introduced as an abbreviation for the square 
bracket in (45), is bounded near the horizon. Hence (58) yields a polynomial 
bound for \fa\. Each multiplication with A~^ and integration increases the 
power by 1 — |, and thus there is a constant ci with 

/ A-^{s)ds j A-^\fa\ < 01^2-^+'' iox p <r < p + e. (61) 

Jo Jo 
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Since 2 — s > 0, (61) is of the order o{u^ ^~^'^), and thus (60) can be written 



as 



a(r) = f A-^s)ds A'^ g + C f A"^ + o{v}-i^'') . (62) 
J p J p J p 

Consider the behavior of the first two summands in (62). The function g 
stands for the last summand in (45). If ii;o 7^ 1, it has a non-zero hmit on the 
horizon. If on the other hand u^o = 1, then g ^ u'^. Substituting into (62) 
and integrating, one sees that the first summand in (62) is ~ u"^ with a given 
by (56). The second summand in (62) vanishes if C = 0, and is ~ u°' with 
a as in (55). According to (53), 1— |<2 — s<2 — s + min(K, s/2). Thus 
the values of a in (55) and (56) are different, and so the first two summands 
in (62) cannot cancel each other. If we choose u so large that 1 — | + > cr, 
(62) yields the Lemma. 

Suppose now that a given by (57) is infinite. Then choosing 



v = max 



1 - |, 2 - s + min(K,s/2) 



we see that the first two summands in (62) are of the order 0{u^) with s 
according to (55) and (56), respectively, and the last summand is of higher 
order. Thus (62) implies that a as defined by (57) is finite (namely, equal 
to the minimum of (55) and (56)), giving a contradiction. Thus a is indeed 
finite. ■ 



In the proof of Proposition 5.4, we already observed that the square 
bracket in the AT^-equation (17) vanishes on the horizon (41). Let us now 
analyze this square bracket in more detail, where we use again the notation 
(51). 

Proposition 7.2. k < 1 and 

[] = Oiu^+'') + Oiu) (63) 

with a as in Lemma 7.1. 



Proof. The derivative of the square bracket is again given by (52). Now 
ao = 0, /3o 7^ and from Lemma 5.1, oq^o + Polo 7^ 0; thus using (21), (22), 
and (54), we get, for r near p, 



[]' = 0(«'^-^+'^) + + 0(1) , 



(64) 
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where we again omitted the expressions "o(u )" and we use the notation 

, . _ J K if 'u;o = 
" \ iiwo^O ■ 

Integrating (64) and using that []r=p = according to (41), we obtain that 
[] = 0(u'^+'^) + + 0{u) . (65) 

Suppose K> 1. Then K + a > 1 and a + (k) + 1 > 1, and (65) becomes 

[] = 0{u). 
We write the ^T^-equation (40) as 

r(^T7 = T3[]-ii^. (66) 

Since {AT'^)' is bounded and = 0{u^^) (by virtue of hypothesis Aj), 
(66) behaves near the event horizon hke 

= 0{u^-^) + 0(n2«-2) . (67) 

Since 2k — 2 > and 1 — ^ < 0, the right side of (67) is unbounded as 
r \ p, giving a contradiction. We conclude that k < 1. 

For K < 1, the second summand in (65) is of higher order than the first, 
and we get (63). ■ 



In the remainder of this section, we shall substitute the power expansions 
(21)-(23) and (54) into the EDYM equations and evaluate the leading terms 
(i.e. the lowest powers in u). This will amount to a rather lengthy consider- 
ation of several cases, each of which has several subcases. We begin with the 
case wq 7^ 0,ibl. The ^-equation (13) simplifies to (46). The ^T^-equation 
(17) for LJ = takes the form (40), and we can for the square bracket use 
the expansion of Proposition 7.2. Finally, we also consider the ^w'-equation 
(16). Using the regularity assumption Aj, we obtain 

yl-eqn: 0{u'-^) = 1 - ~ "^J^^ + 0(u'^) + 0{u'+^''-^) (68) 

e p 

AT^-eqn: n° = 0(u''+"-^) + 0{u^''f) + 0{u^^-^) (69) 

Aw'-eqn: 0{u'+'^-^) = wq {1 - w^) + Oiu'^) + OCu'^^i) . (70) 
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First consider (68). According to An, s + 2k — 2 > 0, and so all powers in (68) 
are positive. We distinguish between the cases where the power s + 2k — 2 
is larger, smaller, or equal to the other powers on the right of (68). Making 
sure in each case that the terms of leading powers may cancel each other, 
we obtain the cases and conditions 

3 

(a) K < S + 2K-2 ^ Wq = 1 ± ep, K = s - 1, s > - (71) 

3 

(b) k = s + 2k-2 =^ wl = l±ep, K = 2- s, s> - (72) 

1 3 

(c) «;>s + 2k;-2>0 =^ wl = 1 ± ep, k = -, s < - (73) 

(d) s + 2k-2 = =^ K = l- |. (74) 

In Case (a), the relations in (71) imply that 

3s 

1-y < 2s - 4 = 2k - 2. 

Hence (69) yields 1 - 3s/2 = ac + cr - 3s/2, so 

cr = 1- K = 2-s . (75) 

This is consistent with Lemma 7.1. But we get a contradiction in (70) as 
follows. Since k = s — 1, we have s + k — 2 = 2s — 3>0;on the other hand, 

a- - = 2-S- - = 2- — < 0. 
2 2 2 

Thus the left-hand side of (70) is bounded, but the right-hand side is un- 
bounded as r\ p. This completes the proof in Case (a) . 

In Case (b), (70) yields that 

^ > I • (76) 

We consider the two cases (55) and (56) in Lemma 7.1. In the first case, 
(76) yields that s < 1, contradicting (53). In the second case, (76) implies 
that s < |. This contradicts the inequality in (72), and thus completes the 
proof in Case (b). 

In Case (c), the relations in (73) give s + k — 2 = s — 1<0, and thus (70) 
implies that s-I-k — 2 = (7 — |, socr=|(s — 1). According to Lemma 7.1, 
(7 = 2 — s or (7 = 1— |. In the first of these cases, we conclude that s = | 
and cr = |. Substituting these powers into (69), we get 



u 



= O(u-i) + 0(u--) + 0{u-^) , 
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which clearly yields a contradiction. Thus cr = 1 — |, giving 

s=-, K = - , cr = - . (77) 

4 2 ' 8 ^ ^ 

This case is ruled out in Lemma 7.3 below. 

In Case (d), wc consider (70). Since s + k — 2 = — k<0, we obtain that 
s + K — 2 = cr— I and thus a = s — 1. Lemma 7.1 yields the two cases 

4 11 

s = -, '^^3' ^^3 ^'^^^ 

3 1 1 

^ = 2' '^=4' ^=2- ^^^^ 

The first of these cases is the exceptional case of Theorem 3.1, and the 
second case is excluded in Lemma 7.4 below. This concludes the proof of 
Theorem 3.1 in the case wq 7^ 0, ±1. 

We next consider the case wp = ±1 . Then the expansions (68)-(70) must 
be modified to 

yl-eqn: 0{u'-^) = 1 + 0{u'^'') + 0{u'+^''-'^) (80) 
^r2-eqn: «° = 0{u^+''-^) + 0{u^-'^) + 0{u'^''-^) (81) 
Aw'-eqn: 0{u'+''-^) = Oiu") + O(tx^-t) . (82) 

One sees immediately that, in order to compensate the constant term in 
(80), s + 2k — 2 must be zero. Hence s + k — 2 = —k<0, and (82) yields 
that s + K — 2 = c7 — I and thus a = s — 1. Now consider Lemma 7.1. In case 
(55), we get the exceptional case of Theorem 3.1, whereas case (56) yields 
that 

3 1 1 

s = - , K = - , a = - . 

2 ' 4 ' 2 

This case is ruled out in Lemma 7.4 below, concluding the proof of Theo- 
rem 3.1 in the case wq = ±1. 

The final case to consider is wq = 0. In this case, the expansions corre- 
sponding to (68)-(70) are 

A-eqn: 0{u'-'^) = 1 - ^ + ©(n^^) + 0{u'+'^''-^) (83) 

e p 

AT'^-eqn: = 0(n'^+'^-^) + Oiu^-^) + 0{u^''-^) (84) 

Aw'-eqn: 0{u'+''-'^) = 0{u'^) + Oiu^'i) . (85) 

If s -h 2av — 2 = 0, we obtain exactly as in the case wq 7^ 0, ±1 above that 
a = s — 1. It follows that k > |, and Lemma 7.1 yields either the exceptional 
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case of Theorem 3.1, or s = 2, contradicting (53). If on the other hand 
s + 2k — 2 > 0, we can in (83) use the inequahty s + 2k — 2 < 2k to conclude 
that s — 1 = s + 2k — 2 and thus k = ^- Now k < §, and Lemma 7.1 together 
with (85) yields the two cases 



s = 
s = 




and 



The first case is ruled out in Lemma 7.3 below, whereas the second case 
leads to a contradiction in (84). This concludes the proof of Theorem 3.1, 
except for the special cases treated in the following two lemmas. 

Lemma 7.3. There is no solution of the EDYM equations satisfying the 
power ansatz (21), (22), and (54) with 

5 13 
^ = 4' = 2' " = 8- 

Proof. Suppose that there is a solution of the EDYM equations with 

A{r) = Aqu^ + o{u^) 
w{r) = wiu^ + o(n2) 

with parameters Ao,wi ^ 0. Consider the A-equation (46). The left side 
is of the order (r — p) 4 . Thus the constant terms on the right side must 
cancel each other. Then the right side is also of the order u^. Comparing 
the coefficients gives 

4^ ° = • 
This equation yields a contradiction because both sides have opposite sign. ■ 



Lemma 7.4. There is no solution of the EDYM equations satisfying the 
power ansatz (21), (22), and (54) with 

3 1 1 

5=2' '^^4' ^ 2' ^^^^ 

Proof. According to (21), we can write the function ^/A as 

y/A = ui + f{u) with / = o(n5). (87) 



F. FINSTER, J. SMOLLER, S.-T. YAU 1255 

Employing the ansatz (21), (22) into the A-equation (46), one sees that 

Aw'"^ = u° + + o{u^) . (88) 
We solve for {w')~^ and substitute (87) to obtain 

1 3 

= ui+u + cif + o{u) (89) 

w 

with a real constant ci. Now consider the ^T^-equation (40), which we write 
again in the form (66) and multiply by A, 

rA(AT^y = (AT^)I ^"^ [] - 4 i^T^) i^w'^) ■ (90) 

As in the proof of Proposition 7.2, a good expansion for the square bracket 
is obtained by integrating its derivative. Namely, according to (52), 

[]' = -w'aP + u'^ + o(uO) 
r 

and thus 

A~^r)[] = - A~^ / {w' a(3){s) ds + + o{u^) . 
f Jp 

Substituting into (90) and using Aj and (88), we obtain 

1 1 1 

A~2{r) / w' aP = + + o{u^) . 
Jp 

We multiply by VA, substitute (87) and differentiate, 

w' a(3 = u"^ + u° + C2 /' + o(n°) . 

Multiplying by 1/w' and using (89) gives the following expansion for a(3, 

aP = + ui + C3 /' + C4U~^ f + o(u^) . (91) 

Next we multiply the ^lu'-equation (16) by VA and write it as 

r'^^/A{VA{VAw')y = e\{VAT){al3) + ui + o(nt) . 

We apply Aj and substitute (87), (88), and (91). This gives an equation of 
the form (modulo higher order terms), 

133 1 133 1 

U2 + -Ui + U4 /' + uf' + U~'i f + = U'^ + Ui + Ui f + U'i f . 
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The constant term ~ v must vanish since ah the other terms tend to zero as 
It — > 0. Furthermore, the terms must cancel because aU the other terms 
are o{u^). We thus obtain 

3 , _1 3 , 

n't f + u 4 / = u4 -\- u f + f , 

so that 

u f' + f = u + ui f' + f = u + o{u) , 

and we find that / satisfies an equation of the form 

di u f' + d2 f = d^u + o{u) . 

A straightforward but tedious calculation yields that the coefficients di and 
^2 both vanish, and that ds is non-zero. This is a contradiction. ■ 



8 The Exceptional Case 



In this section, we consider the exceptional case 

4 1 1 

3 ' 3 ' 3 

By employing the power ansatz (21), (22), and (54) into the EDYM equa- 
tions and comparing coefficients (using Mathematica) , we find that the 
solution near the event horizon is determined by the five free parameters 
(/3o, 70, "1, c, p). The remaining parameters are given by 



Ao ^ 




•^o = \h'^-Pi + — Poio 

" rm 



Wo 



9Po \rin?r'^ Pi - 2cmr Polo + 7o 



r Y — (1 — Wq)'^ 

mrPo - c7o 



Wl 



rVAo 
9 aoPo 
2r Ao 



Expanding to higher order, we obtain after an arduous calculation two fur- 
ther constraints on the free parameters, thus reducing the problem to one 
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involving only three parameters. We investigated this three-par amctcr space 
numerically starting at vq = p + e and found strong evidence that no global 
black hole solutions exist. Indeed, either the power ansatz was inconsistent 
near the event horizon (that is, for r close to ro + £ the numerical solution 
deviated from the power ansatz, and became singular as £ \ 0), or else the 
solutions developed a singularity for finite r. 
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